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Abstract 

The theory of derivatives and integrals of non-integer order goes back to Leibniz, Liou- 
ville, Grunwald, Letnikov and Riemann. The history of fractional vector calculus (FVC) 
has only 10 years. The main approaches to formulate a FVC, which are used in the physics 
during the past few years, will be briefly described in this paper. We solve some problems 
of consistent formulations of FVC by using a fractional generalization of the Fundamen- 
tal Theorem of Calculus. We define the differential and integral vector operations. The 
fractional Green's, Stokes' and Gauss's theorems are formulated. The proofs of these the- 
orems are realized for simplest regions. A fractional generalization of exterior differential 
calculus of differential forms is discussed. Fractional nonlocal Maxwell's equations and 
the corresponding fractional wave equations are considered. 
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1 Introduction 



The fractional calculus has a long history from 30 September 1695, when the derivative of 
order a = 1/2 has been described by Leibniz [H[2] (see also [6]). The theory of derivatives and 
integrals of non-integer order goes back to Leibniz, Liouville, Grunwald, Letnikov and Riemann. 
There are many interesting books about fractional calculus and fractional differential equations 
[H El El HI [5] (see also [8]). Derivatives and integrals of fractional order, and fractional 
integro-differential equations have found many applications in recent studies in physics (for 
example, see books P [TDJ QU Q2J, and reviews [El (H [15]). 

The history of fractional vector calculus (FVC) is not so long. It has only 10 years and 
can be reduced to the papers PS]- [27]. The main approaches to formulate a FVC, which are 
used in the physics during the past few years, will be briefly described in this paper. There are 
some fundamental problems of consistent formulations of FVC that can be solved by using a 
fractional generalization of the Fundamental Theorem of Calculus. Fractional vector calculus 
is very important to describe processes in fractal media (see for example [10J). A consistent 
FVC can be used in fractional electrodynamics [HI HHJ [20], [25] and fractional hydrodynamics 

[mug. 

In Section 2, we describe different approaches to formulate FVC, which are used in the 
physics during the past 10 years. The problems of consistent formulation of FVC are described 
in Section 3. A fractional generalization of the Fundamental Theorem of Calculus is considered 
in Section 4. In Section 5, the differential and integral vector operations are defined. In Sections 
6-8, the fractional Green's, Stokes' and Gauss's theorems are formulated. The proofs of these 
theorems are realized for simplest regions. In Section 9, a fractional generalization of exterior 
calculus of differential forms is discussed. In Section 10, fractional nonlocal Maxwell's equations 
and the corresponding fractional wave equations are considered. 
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2 Approaches to fractional vector calculus 



For Cartesian coordinates, fractional generalizations of the divergence or gradient operators 
can be denned by 

grad a f{x) = e s D a J\x) 1 (1) 

div a F(x) = D°F s (x), (2) 

where Df are fractional (Liouville, Riemann-Liouville, Caputo, etc.) derivatives [H [21 El IH [5] 
of order a with respect to x s , (s = 1,2,3). Here e s (s = 1,2,3) are orthogonal unit vectors, 
and F s (x) are components of the vector field 

F(x) = F s (x)e s = F x e x + F y e y + F z e z . (3) 

The main problem of formulation of FVC appears, when we try to generalize the curl 
operator and the integral theorems. In Cartesian coordinates, the usual (integer) curl operator 
for the vector field ([3]) is defined by 

cur IF = eiei mn D m F n , (4) 

where D m = d/dx m , and E\ mn is Levi-Civita symbol, which is 1 if k) is an even permutation 
of (1,2,3), (—1) if it is an odd permutation, and if any index is repeated. The Fourier 
transform of the curl operator is 

T(curlF(x)) = eiei mn (ik m )F n (k), (5) 

where 

/+oo 
<('•■>■< lh 'F„(r). (6) 
-co 

To define a generalization of PJ, we can use a fractional integro-differentiation instead of the 
derivative D m . 
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2.1 Ben Adda's fractional vector calculus 

In the paper [IB] (see also [17]). fractional generalizations of gradient, divergence and curl 
operator for analytic functions have been suggested in the form 

grad a f(x) = \ e s D«f(x), (7) 
1 {a + 1) 



1 (a + 1) 



curl a F = - 1 —- ei e lmn D^F n (x), (9) 
1 [a + 1) 

where T(a + 1) is the Gamma function. In these definitions, the Nishimoto fractional derivative 
[8] (see also Section 22 of [2]) is used. This derivative is a generalization of the Cauchy's 
differentiation formula. 

Fractional generalizations of integral operations (flux and circulation), and generalizations 
of Gauss's, Stokes', Green's integral theorems are not considered. 

2.2 Engheta's fractional vector calculus 

In the paper [18] (see also [HI [191 120]), a fractional generalization of curl operator has been 
suggested in the form 

curl a F = e l e lrnn D°F n (x), (10) 

where are fractional Liouville derivatives [5] of order a with respect to x m , (m = 1,2,3), 
that are defined by 

D a m f{x):= lim a D a x J(x m ). (11) 

a — >— oo 

Here a D x is the Riemann-Liouville derivative 

i o n r f(x') 

a D a J(x) = -— / J -±-^—dx', (n-Ka<n). (12) 

xJ K ' T{n- a)dx n J a (x - x') a ~ n+1 

The fractional Liouville derivative (TTTj) can be defined through the Fourier transform by 

D a m F n {x)=T- x [{ik m ) a F n {k)) = —--J d 3 ke ikx (tk m ) a F n (k), (13) 



where F n (k) is defined by (EI) and i a = exp{za7r sgn(k) /2} . For this fractional curl operator, 
the fractional integral Stokes' and Green's theorems are not suggested. The problems of a 
generalization of these theorems will be considered in the next section. 

In general, the fractional vector calculus must include generalizations of the differential oper- 
ations (gradient, divergence, curl), the integral operations (flux, circulation), and the theorems 
of Gauss, Stokes and Green. 

2.3 Meerschaert-Mortensen-Wheatcraft fractional vector calculus 

In the paper of Meerschaert, Mortensen and Wheatcraft [21], a fractional generalization of curl 
operator has been suggested as 

curPF = ^e lmn D r Jl- a F n , (14) 

where I^~ a are fractional integrals of order (1 — a) with respect to x n , (n = 1, 2, 3). Note that 
the integration I\~ a in (1141) is considered with the index n as the component F n . The derivative 
D m = d/dx m in Eq. ffl~4"l) is considered with respect to x m , where m ^ n. Therefore expression 
(}T4"|) can be presented as the usual (integer) curl operator 

curl a F = curlF(a) (15) 

for the field 

F(a) = eJl^Fn. (16) 

Equation ffl5|) allows us to use the usual (integer) integral Stokes' and Green's theorems. 

In Eq. ( fl4l . the fractional integral 1° and the integer derivative D m have antisymmetric 
indices, and the components of (fl4l) are 

{currF) x = D y l\~ a F z - D z F y - a F y , (17) 

(curl a F) y = D z l\~ a F x - DJl- a F Z} (18) 
(curl a F) z = DJ l y - a F y - D y F- a F x . (19) 
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It is easy to see that operator (|T4l) has no fractional derivatives with respect to x m , (m — 1, 2, 3), 
like as D a m = D m Ii~ a or C D^ = I^ a D m . 

As a result, we have the usual (integer) vector calculus for new type of fields as in (|T6|) . The 
suggested approach cannot be considered as a fractional generalization of vector calculus. It is 
important to define a curl operator with fractional derivatives in such a form that fractional 
generalizations of the integral theorems exist. 

2.4 Other approaches to fractional vector calculus 

In the papers [251 121], fractional generalizations of integral operations and Gauss's, Stokes', 
Green's theorems have been suggested. These generalizations are considered to describe frac- 
tional media by a continuous medium model. The differential operations are defined with 
respect to fractional powers of coordinates. These operations are connected with fractional 
derivatives only by Fourier transforms (see [28J ) . As a result, an "ideal" fractional vector 
calculus is not suggested. 

In the papers [261 12Z] , fractional differential vector operations are considered by using frac- 
tional generalizations of differential forms that are suggested in [22] (see also [231 [23]). A frac- 
tional gradient is defined by an exact fractional 1-form. A fractional curl operator is described 
by a fractional exterior derivative of a fractional differential 1-form. The Riemann-Liouville 
derivatives are used in [26J, and the fractional Caputo derivatives are used in [27]. We have 

grad a f(x)=e s ^Dlf(x), (20) 

curl a ¥ = e l e lmn CD a Xm F n} (21) 
where oD% m is a fractional Caputo derivative with respect to x m \ 



c a D«f(x) = — / —dx' J \ ' , (n - 1< a < n). (22) 

a xJy ' Y(n-a)J a (x-x') a - n+1 d(x') n v ) \ ) 

The fractional generalizations of integral theorems (Gauss's, Stokes', Green's theorems) are not 
considered and the fractional integrals for differential forms are not defined. 
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3 Problems of fractional generalization of Green's for- 
mula 



Let us describe a main problem that appears when the curl operator and integral formulas are 
generalized on a fractional case. For simplification, we consider a rectangular domain on IR 2 
and integral formulas in Cartesian coordinates. 
The Green's formula in Cartesian coordinates is 



(23) 



/ (F x dx + F y dy)= [ [ dxdy [D y F x — D x F y ] , 

JdW J JW 

where F x = F x (x,y) and F y = F y (x,y) are functions defined for all (x, y) in the region W. 
Let W be the rectangular domain 

W := {(x, y) : a < x <b, c < y < d} 

with the sides AB, BC, CD, DA, where the points A, B, C , D have coordinates 

A(a,c), B(a,d), C{b,d), D(b,c). 

These sides form a boundary dW of W . Then 

(F x dx + Fydy) = / F x dx + / F x dx + / F y dy + / F y dy = 
J bc Jda Jab Jcd 

pb pa pd pc 

i F x (x,d)dx+ / F x (x, c)dx+ / F y (a,y)dy+ / F y {b,y)dy = 

a J b J c J d 

rb fd 

- dx[F x (x,d)-F x (x,c)}+ dy[F y (a,y)-F y (b,y)]. (24) 

J a J c 



The main step of proof of Green's formula is to use the Newton-Leibniz formula 



dxDJ(x) = f(b)-f(a). 



(25) 



The function f(x) in (125]) is absolutely continuous on [a, b}. As a result, expression fT2"4"|) can be 
presented as 



dx 



dy D y F x (x,y) 



+ J dy 
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dx D x F y (x,y) 



= dx dy [D y F x (x,y) - D x F y (x,y)} = / / dxdy [D y F x — D x F y ] . 

Ja Jc J JW 

To derive a fractional generalization of Green's formula (I23I) . we should have a generalization 
of the Newton-Leibniz formula (125]) in the form 

X a D:f(x) = f(b)-f(a), (26) 

where some integral and derivative of noninteger order are used. This generalization exists for 
specified fractional integrals and derivatives, and does not exist for arbitrary taken type of the 
fractional derivatives. 

For the left Riemann-Liouville fractional integral and derivative (Lemma 2.5. of [5]), we 
have 

a i? a D a x f(x) = f(b) - r f ~ Q) T' pr j a/rvxa), (27) 

where D x ~^ = d n ~i / 'dx n ~i are integer derivatives, and n— 1 < a < n. In particular, if < a < 1, 
then 

JZaD«J{x) = fib) - {b ~ fl)Q " 1 Jl- a f(x), (28) 

r{a) 

Obviously that Eq. (1281) cannot be considered as a realization of (1261) . The left Riemann- 
Liouville fractional integral for x e [a, b] is defined by 



1 r dr' 

J a J{x) := — - / aX ,. n _ a (a > 0). (29) 
r (a) L (x — x ni a 



The left Riemann-Liouville fractional derivative for x G [a, b] and n — l<a<nis defined by 
Note that Eq. (|27|) is satisfied if /(x) is Lebesgue measurable functions on [a, b] for which 



/ /(x) <ix < oo, 

■/ a 



and a I^ a f(x) of the right-hand side of ( 1271) has absolutely continuous derivatives up to order 
(n — 1) on [a, b}. 
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Properties (}2Tj) and ( 1281) are connected with the definition of the Riemann-Liouville frac- 
tional derivative, where the integer-order derivative acts on the fractional integral: 

a D a x = D n x a /r°, (n - 1< a < n). (31) 

This definition gives that the left-hand side of (1281) is 

ja r~)d jet 7~jn jn—ct ( e \ r )\ 

a L x a u x — a L x u x a L x i 

where the integer derivative D™ is located between the fractional integrals. Since the operations 
D™ and a J° are not commutative 

a L x u x u x a L x T U ) 

we get the additional terms, which cannot give the right-hand side of (I26I) . This noncommu- 
tativity can be presented as a nonequivalence of Riemann-Liouville and Caputo derivatives 

BE], 

rl_ 1 / \j—Q, 

C a D a J{x) = a D a J(x) - r )~ a *^ (Dif)(a), (n-Ka<n). (33) 

The left Caputo fractional derivative is defined by the equation (compare with fl3~Tj) ) 

c a D«f{x) := j:~ a D n J(x), (n - K a < n). (34) 

The noncommutativity of D x and a I x in (I3"2I) does not allow us to use semi-group property (see 
Lemma 2.3 of [5] and Theorem 2.5 of [2]) of fractional integrals 

a ^ X = a i: +f3 , (a > 0, (3 > 0). (35) 

Note that equation (1351) is satisfied at almost every point x G [a, 6] for f(x) G L p (a,b) and 
a, /3 > 0. We denote by L p (a,b) (1 < p < oo) the set of those Lebesgue measurable functions 
on [a, b] for which 

J dx\f(x)\ p j <oo. 



In general, the semi-group property 

a D a x a Dl = a D^, (a > 0, 0). (36) 
9 



is not satisfied for fractional derivatives (see Property 2.4 in [5]). For some special cases, Eq. 
fl36|) can be used (see Theorem 2.5. in [2J). For example, the property fl36|) is satisfied for the 
functions 

f(x)e J^{L x {a,b% 
i.e., equation ( 1361) is valid for f(x) if there exists a function g(x) G L\(a,b) such that 

/(*) = a i: +p g(x). 

The semi-group property for fractional derivatives is also valid if a = 0, b = oo and /(x) is 
infinitely differentiable (generalized) function on [0, oo) (see Sec. 1.4. 5. of [30] and Sec. 8. 3. of 
12]). 

In order to have a fractional generalization of the Newton-Leibniz formula of the form ( 1261) . 
we must replace the left Riemann-Liouville derivative a D£ in Eq. (!26l) . where 

ja 7-ja ja/-r\n jn—o.\ 

by the left Caputo derivative „ D", such that the left-hand side of (|26|) is 
Then, we can use the semi-group property (1351) . and 

a /: c a D a j{ X ) = j: a ir a D n j(x) = „/;.7r/ ./•(,■). 

In particular, if n = 1 and < a < 1, then 

a / b Q ?££/(*) = ^^/(x) = fdxD l J{x) = f(b) - f(a). (37) 

As a result, to generalize Gauss's, Green's and Stokes' formulas for fractional case, we can 
use the equation with the Riemann-Liouville integral and the Caputo derivative: 

*I?Zr%f{x) = f{b)-f(a). (38) 
This equation can be considered as a fractional analog of the Newton-Leibniz formula. 
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4 Fractional Generalization of the Fundamental Theo- 
rem of Calculus 



The fundamental theorem of calculus (FTC) is the statement that the two central operations 
of calculus, differentiation and integration, are inverse operations: if a continuous function is 
first integrated and then differentiated, the original function is retrieved 

Dl Jlf{x) = /(*). (39) 

An important consequence, sometimes called the second fundamental theorem of calculus, al- 
lows one to compute integrals by using an antiderivative of the function to be integrated: 

J 1 b D 1 J(x) = f(b)-f(a). (40) 

If we use the Riemann-Liouville integrals and derivatives [21 [5] , we cannot generalize (j4"0~j) for 
fractional case, since 

Jt aDZf(x)?f(b)-f(a), (41) 

In this case, we have equation (1271) . 

The FTC states that the integral of a function / over the interval [a, b] can be calculated by 
finding an antiderivative F, i.e., a function, whose derivative is /. Integral theorems of vector 
calculus (Stokes', Green's, Gauss's theorems) can be considered as generalizations of FTC. 

The fractional generalization of the FTC for finite interval [a, b] can be realized (see remarks 
after proof of the theorem and Section 3.) in the following special form. 

Fundamental Theorem of Fractional Calculus 

(1) Let f{x) be a real-valued function defined on a closed interval [a, b] . Let F(x) be the function 
defined for x in [a, b] by 

F(x) = j:f(x), (42) 
where a I£ is the fractional Riemann-Liouville integral 

^■■=W)S1j^w=-^ (43) 
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then 

c a D a x F{x) = f(x) (44) 
for x G (a, b), where is the Caputo fractional derivative 

1 r x A T > r\ n F{r\ 

° a D«F{x) = J n x - a D n x F(x) = / -, ^— (n - 1< a < n). (45) 

a x w x x v ' T(n- a) J a (x - x') 1+a ~ n dx n ' v ' 

(2) Let f(x) be a real-valued function defined on a closed interval [a,b]. Let F{x) be a 

function such that 

f(x) = c a D a x F{x) (46) 

for all x in [a, b], then 

JZf(x) = F(b)-F(a), (47) 

or, equivalently, 

° a D«F{x) = F{b) - F{a), (0 < a < 1). (48) 

As a result, we have the fractional analogs of equations (|39l) and (|4*0|) in the form 

c a D a x J£f(x) = f(x), (a > 0), (49) 

a ^ C a D«F{x) = F(x) - F(a), (0 < a < 1), (50) 
where a I x is the Riemann-Liouville integral, and ^D x is the Caputo derivative. 

Proof. 

The proof of this theorem can be realized by using the Lemma 2.21 and Lemma 2.22 of [5]. 

(1) For real values of a > 0, the Caputo fractional derivative provides operation inverse to 
the Riemann-Liouville integration from the left (see Lemma 2.21 [5]), 

c a D« a I«f{x) = f{x), (a>0) (51) 

for f(x) G L^a^b) or f(x) G C[a,b]. 

(2) If f(x) G AC n [a,b] or f(x) G C n [a,b], then (see Lemma 2.22 [5]) 

n— 1 

a ^ C a D a x f{x) = f(x) -£-(*- aY(l)>J)(a). (n-l<a<n), (52) 



3=0 J 
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where C n [a, b] is a space of functions, which are n times continuously differentiable on [a, b]. In 
particular, if < a < 1 and f(x) G AC[a, b] or f(x) G C[a, b], then 

J« c a D«f{x)=f{x)-f{a). (53) 

This equation can be considered as a fractional generalization of the Newton-Leibniz formula 
in the form f[2"B"l) . 

Remark 1. In this theorem (see Eqs. ( 1421148]) ). the spaces Li[a, b] and AC[a, b] are used. 

(a) Here AC [a, b] is a space of functions F(x), which are absolutely continuous on [a, b]. It 
is known that AC [a, b] coincides with the space of primitives of Lebesgue summable functions 
and therefore an absolutely continuous function F(x) has a summable derivative D].{x) almost 
everywhere on [a, b]. If F(x) G AC[a, b], then the Caputo derivative (0 < a < 1) exists almost 
everywhere on [a, b] (see Theorem 2.1 of [5]). 

(b) We denote L p (a, b) the set of those Lebesgue measurable functions / on [a, b] for which 

\\f\\p=(j'\m\ p dx^ /P <oo. (54) 

If f(x) G L p (a, b), where p > 1, then the fractional Riemann-Liouville integrations are bounded 
in L p (a,b), and the semi-group property 

a ^ jPf(x) = a i: + ?f(x), (a > 0, (5 > 0) (55) 

are satisfied at almost every point x G [a, b\. If a + (3 > 1, then relations (|55|) holds at any 
point of [a, 6] (see Lemma 2.1 and Lemma 2.3 in [5]). 

Remark 2. For the Riemann-Liouville derivative a D", the relation 

a D a x JZfix) = /(z), (a > 0) (56) 

holds almost everywhere on [a, 6] for /(x) G L p (a,b) (see Lemma 2.4 of [3]). 

Remark 3. The Fundamental Theorem of Fractional Calculus (FTFC) uses the Riemann- 
Liouville integration and the Caputo differentiation. The main property is that the Caputo 
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fractional derivative provides us an operation inverse to the Riemann-Liouville fractional in- 
tegration from the left. It should be noted that consistent fractional generalizations of the 
FTC, the differential vector operations and the integral theorems for other fractional integro- 
differentiation such as Riesz, Grunvald-Letnikov, Weyl, Nishimoto are open problems. 



Remark 4. In the theorem, we use < a < 1. As a result, we obtain the fractional 



Green's, Stokes' and Gauss's theorems for < a < 1. Equation (f49l) is satisfied for a G R+. 
The Newton-Leibniz formula fl50l) holds for < a < 1. For a > 1, we have fl52l) . As a result, 
to generalize the Green's, Stokes' and Gauss's theorems for a G R + , we can use Eq. ( )52|) in the 
form 

n— 1 

/(&) - /(a) = a / b Q c a D a J(x) + £)-(&- «)V (i) (a)> (n - K a < n), (57) 

i=i J ' 

where f^\x) = D J x f(x). In particular, if 1 < a < 2, then n = 2 and 

/(&) - f(a) = Jt C a D a J{x) + (b - a)V'(a). (58) 

Remark 5. In the FTFC, we use the left fractional integrals and derivatives. The Newton- 
Leibniz formulas can be presented for the right fractional Riemann-Liouville integrals and the 
right fractional Caputo derivatives in the form 

x c x D%f{x) = /(*) - \ ( b - *y- ( 59 ) 

In particular, if < a < 1, then 

J^D«f(x) = f(x)-f(b). (60) 

For Im(a) = 0, a > 0, /(x) G -^^(a, 6) or /(x) G C[a, b], then 

?D?«if/(x) = /(x). (61) 

As a result, fractional generalization of differential operations and integral theorems can be 
defined for the right integrals and derivatives as well as for the left ones. 
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5 Definition of fractional vector operations 
5.1 Fractional operators 

To define fractional vector operations, we introduce the operators that correspond to the frac- 
tional derivatives and integrals. 

We define the fractional integral operator 

J >'^W)Lv^' (a>0) ' (62) 

which acts on a real-valued function f(x) G Li[a, b] by 

^y^=w)Slv^f-- (63) 

The Caputo fractional differential operator on [a, b] can be defined by 

c a D a x [x'] := / (n - 1< a < n), (64) 

a xL 1 T(n-a)J a (x - x') 1+a - n dx' n v ' v 1 

such that the Caputo derivatives for f(x) G AC n [a, b] is written as 



c a D a \x']f(x') = , r %- ^A^, (n-l<a<n). (65) 

" 1 1 T{n-a)J a (x-x') l+a - n dx' n V 



It is easy to see that 

c a D a x [x'\ = a r~ a \x'\ D [x'\ , (n - 1< a < n) . 
Using these notations, formulas ( H9l) and ([50ll of the FTFC can be presented as 

C a D a x {AJ a x ,[x"]f{x") = f{x\ («>0), (66) 

J«[x] ^D a x [x']f(x') = f{b) - f(a), (0<a< 1). (67) 

This form is more convenient than f|4*9|) and fl50|) . since it allows us to take into account the 
variables of integration and the domain of the operators. 
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5.2 Definition of fractional differential vector operations 

Let us define the fractional differential operators (grad, div, curl) such that fractional gener- 
alizations of integral theorems (Green's, Stokes', Gauss') can be realized. We use the Caputo 
derivatives to defined these operators and we use the Riemann-Liouville integrals in the gener- 
alizations of the integral theorems. 

Let W be a domain of M 3 . Let f(x) and F(x) be real- valued functions that have continuous 
derivatives up to order (n — 1) on W, such that the [n — 1) derivatives are absolutely continuous, 
i.e., /, F e AC"[W]. We can define a fractional generalization of nabla operator by 

= C D^ = ei C D^[x\ + e 2 c D^[y\ + e 3 c D&[z], (n - 1< a < n). (68) 

Here, we use the fractional Caputo derivatives c Dy V [x m ] with respect to coordinates x m . For 
the parallelepiped 

W := {a < x < b, c < y < d, g < z < h}, 

we have 

c Dfr[x]=?L%[x], C D^[y}= c c D«[yl c ' D%[z] = c g D« h [z]. (69) 
The right-hand sides of these equations the Caputo derivatives are used. 

(1) If / = f(x, y, z) is (n — 1) times continuously differentiate scalar field such that -D™" 1 / 
is absolutely continuous, then we define its fractional gradient as the following 

Grad^f = c B^f = e, c D^[ Xl ]f(x, y, z) = 

= ei c D^[x]f(x,y,z) + e 2 c D^[y]f(x,y, z) + e 3 c D^[z]f(x,y, z). (70) 

(2) If F(x, y, z) is (n — 1) times continuously difierentiable vector field such that D™~ l Fi are 
absolutely continuous, then we define its fractional divergence as a value of the expression 

Div^F = ( C D^,F) = c D^[xi\Fi(x,y,z) = 

= c D^[x]F x (x,y,z)+ c D^[y]F y (x,y,z)+ c D^[z]F z (x,y, z). (71) 
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(3) The fractional curl operator is defined by 



Curl^F 



Gie lnk c D^[x m ]F k = ex ( C D^[y]F z - c D^[z\F y ) + 



+ e 2 ( c D^[z}F x - c D^[x]F z ) + e 3 ( C D%[x]F v - c D^[y]F x ) , (72) 

where F k = F k (x, y, z) e AC n [W], (k = 1, 2, 3). 

Note that these fractional differential operators are nonlocal. As a result, the fractional 
gradient, divergence and curl depend on the region W. 

5.3 Relations for fractional differential vet or operations 

(a) The first relation for the scalar field / = f(x, y, z) is 

Curl^ Grad a w f = e t e lmn c D^[x m ] c Bfr[x n )f = 0, (73) 

where Ei mn is Levi-Civita symbol, i.e. it is 1 if (i,j,k) is an even permutation of (1,2,3), (— 1) 
if it is an odd permutation, and if any index is repeated, 
(b) The second relation, 

3 

Div^Grad^f(x,y,z) = c D^[ Xl ] c 'D^fix^z) = ^2( c D^[ Xl ]) 2 f(x,y,z). (74) 



i=i 



Using notation fl68l) . 



In the general case, 



It is obvious from 



Divyy Grad^ = ( C D^) 2 = ( C D^, C D^) . (75) 
( C Z^N)V C D 2 M. (76) 



\a U x! ~ a L x U x a 1 x U x 

~~ a 1 x a L x U x U x ' a L x \ U x 1 a L x \ U x ~ a U x ~>~ a L x \ u x i a L x \ U . 

where 

[nil rn-ai r^n jn — a jn — a r-jra 7-jQ C T~\a I r\ 
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(c) It is easy to prove the following relation, 

Divw Curl^F(x, y, z) = c D^[xi}ei mn c D^[x m ]F n (x, y, z) = 

= e lmn c D^[ Xl ] c D^[x m }F n {x, y, z) = 0, (77) 

where we use antisymmetry of ei mn with respect to m and n. 

(d) There exists a relation for the double curl operation in the form 

Curlw CurlyyF(x, y, z) = eieimn C D^ v [x m ]e npq c D^[x p ]F q (x,y, z) = 

= ^l^lmn^npq Dy^\Xm\ D^r[Xp]Fq(x, y , z) . (^8) 

Using 

E-lmn^-lpq &mp$nq &mq&npi (79) 

we get 

Curl^ Curl^ F(x, y, z) = Grad a w Div^ F(x, y,z)-( c B^) 2 F(x, y, z), (80) 

(e) In the general case, 

c a mA(mg{x')) (si) 

For example (see Theorem 15.1. from [2]), 

*m\A (/OOgOO) = fj T (j + ijr(a -j + 1) ( aD "~ j W f ^) ( D * 9 ^) > (82) 

if f(x) and g(x) are analytic functions on [a, b\. As a result, we have 

Grad a w (/«?) ± (Grad a w f) g + (Orad^ f, (83) 

Divfr (/f) ^ (Grad^f, f) + / DivfrF. (84) 

These relations state that we cannot use the Leibniz rule in a fractional generalization of the 
vector calculus. 
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5.4 Fractional integral vector operations 

In this section, we define fractional generalizations of circulation, flux and volume integral. 
Let F = F(x, y, z) be a vector field such that 

F(x, y, z) = eiF x (a;, y, z) + e 2 F y (x, y, z) + e 3 F z (x, y, z). 

If F X) F y , F x are absolutely integrable real-valued functions on M 3 , i.e., F X) F y , F x e Li(R 3 ), 
then we can define the following fractional integral vector operations of order a > 0. 

(1) A fractional circulation is a fractional line integral along a line L that is defined by 

S2(F) = (l% F) = I*[x\F x + I2[y}F y + I a L [z}F z . (85) 

For a — 1, we get 

4(F) = (li,F) = ^(dL,F) = J(F x dx + F v dy + F g dz), (86) 
where dL = e^rr + e 2 dy + e 3 rfz. 

(2) A fractional flux of the vector field F across a surface S is a fractional surface integral 
of the field, such that 

S§(F) = (l£,F) =/|[ J /,^ se + /^,x]F w + /|[x,i/]iJ;. (87) 

For a = 1, we get 

$i(F) = (ls,F) = (<*S,f) = J J (F x dydz + F y dzdx + F z dxdy), (88) 
where dS = e\dydz + e 2 dzdx + e 3 dxdy. 

(3) A fractional volume integral is a triple fractional integral within a region in IR 3 of a 
scalar field / = f(x, y, z), 

W) = /^[x, y, y, *) = 7^[x]7^[y]7^[z]/(x, y, z). (89) 

For ct = 1, we have 

nU/):= / // dVf(x,y,z)= Iff dxdydz f(x,y, z). (90) 
J J Jw J J Jw 

This is the usual volume integral for the function f(x,y,z). 
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6 Fractional Green's formula 



Green's theorem gives the relationship between a line integral around a simple closed curve dW 
and a double integral over the plane region W bounded by dW. The theorem statement is the 
following. Let dW be a positively oriented, piecewise smooth, simple closed curve in the plane 
and let W be a region bounded by dW. If F x and F y have continuous partial derivatives on an 
open region containing W , then 

J (F x dx + F y dy}= J J (p y F x -D x F^)dxdy. (91) 

A fractional generalization of the Green's formula ( |9T|) is presented by the following statement. 

Theorem (Fractional Green's Theorem for a Rectangle) 
Let F x (x,y) and F y (x,y) be absolutely continuous (or continuously differentiable) real-valued 
functions in a domain that includes the rectangle 

W := {{x, y) : a < x < b, c<y <d}. (92) 

Let the boundary ofW be the closed curve dW . Then 

I^ w [x}F x (x,y) + I^ w [y}F y (x,y) = I^[x,y] ( c D$ w [y]F x (x,y) - c D% w [x]F y (x,y)) , (93) 

where < a < 1. 

Proof 

To prove equation (1931 . we change the double fractional integral /^[x, y] to the repeated frac- 
tional integrals I^[x] Iw[y], and then employ the Fundamental Theorem of Fractional Calculus. 

Let W be the rectangular domain (|92|) with the sides AB, BC, CD, DA, where the points 
A, B, C, D have coordinates 

A(a,c), B(a,d), C{b,d), D(b,c). 

These sides form the boundary dW of W. 
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For the rectangular region W defined by a < x < b, c < y < d, the repeated integral is 

Iw\A Iwiy] = al^[x} Jd [y] > 

and equation ( |93|) is 

JZ{x} (F x (x,d) - F x (x,c)) + J2[y] (F y (a,y) - F y (b,y)) = 

= am J a d [y] ( °D«[y>)F x (x,y') - ° D«[x')F y (x> , y)) . (94) 
To prove of the fractional Green's formula, we realize the following transformations 

(l a aw , F) = I% w [x]F x + I% w [y]F y = r BC [x]F x + I% A [x]F x + I% B [y}F y + I^ D [y]F y = 

= a mF x (x,d) - a I?[x]F x {x,c) + c I^[y]F y (a,y)dy- c I^[y}F y (b,y) = 

= a m [Fx(x,d) - F x (x,c)} + J2\V) [F y (a,y) - F y {b,y)). (95) 
The main step of the proof of Green's formula is to use the fractional Newton-Leibniz formula 

F x {x, d) - F x (x, c) = c I«[y] c c D«[y'\ F(x, y% 

F y (a, y) - F y (b, y) = - c a D«[x'\ F(x', y). (96) 

As a result, expression fl95|) can be presented as 

a m {jM C c D«[y'\F x {x,y')) + c I a d [y] {- J?[x] ° a D a x [x'] F y (x', y)} = 

= am cm { C cD a y [y'\F x {x,y') - c a D«[x')F y {x\y)) = 
= Ifr[x,y] ^D^y'}F x (x,y') - c a D«[x'] F y (x', y)) . 
This is the left-hand side of Eq. (1941) . This ends the proof. 

Remark 1. In this fractional Green's theorem, we use the rectangular region W. If the 
region can be approximated by a set of rectangles, the fractional Green's formula can also be 
proved. In this case, the boundary dW is presented by paths each consisting of horizontal and 
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vertical line segments, lying in W. 



Remark 2. To define the double integral and the theorem for nonrectangular regions R, we 
can consider the function f(x, y), that is defined in the rectangular region W such that R C W 
and 

{F(x,y), (x,y) e R; 
(97) 
0, (x,y)eW/R. 

As a result, we define a fractional double integral over the nonrectangular region R, through 

the fractional double integral over the rectangular region W: 

I%[x,y]F(x,y)=T$r[x,y]f(x,y). (98) 

To define double integrals over nonrectangular regions, we can use a fairly general method 
to calculate them. For example, we can do this for special regions called elementary regions. 
Let R be a set of all points (x, y) such that 

a < x < b, ipi(x) < y < ^{x). 

Then, the double integrals for such regions can be calculated by 

l a R [x,y} F(x,y) = a l a b [x] Vl(x) l« 2(x) [y} F(x,y). (99) 

It is easy to consider the following examples. 

1) (pt(x) = 0, y = cp 2 (x) = x 2 , F(x, y) = x + y. 

2) ip\(x) = x 3 , (fi(x) = x 2 , F(x, y) = x + y. 

3) ifi(x) = 0, y = (p 2 (x) = x, F(x,y) = xy. 

The fractional integrals can be calculated by using the relations 

a i:[x}(x - af = - Jg±iL ( g - af +a , (100) 
1 [p + a + 1) 

where a > 0, (3 > 0. For other relations see Table 9.1 in [2]. To calculate the Caputo derivatives, 
we can use this table and the equation 

° a D<* x [x']f{x>) = a D«[x>]f{x') - £ r / fc)(a | lV n - 1< a < n. (101) 
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Note that the Mittag-Leffler function E a [(x' — a) a ] is not changed by the Caputo derivative 

c a D a x [x'} E a [(x' - a) a ] = E a [(x - a) a }. (102) 

This equation is a fractional analog of the well-known property of exponential function of the 
form D\ exp(x — a) = exp(x — a). Therefore the Mittag-Leffler function can be considered as 
a fractional analog of exponential function. 



7 Fractional Stokes' formula 

We shall restrict ourselves to the consideration of a simple surface. If we denote the boundary 
of the simple surface W by dW and if F is a smooth vector field defined on W, then the Stokes' 
theorem asserts that 



F,dL) = / (curlF,dS). (103) 

The right-hand side of this equation is the surface integral of curlF over W, whereas the left- 
hand side is the line integral of F over the line dW. Thus the Stokes' theorem is the assertion 
that the line integral of a vector field over the boundary of the surface W is the same as the 
integral over the surface of the curl of F. 
For Cartesian coordinates, Eq. (11031) gives 



dW 



(104) 



(F x dx + Fydy + F z dz) = 

■ \ / 

= J J (dydz [D y F z - D z F y ] + dzdx [D Z F X - D X F Z ] + dxdy [D x F y - D y F x 
Let F = F(x, y, z) be a vector field such that 

F(x, y, z) = etF x (x, y, z) + e 2 F y (x, y, z) + e 3 F z (x, y, z). 

where F x , F y , F x are absolutely continuous (or continuously differentiable) real- valued functions 
on R 3 . Then the fractional generalization of the Stokes' formula (11041) can be written as 

(l a dw ,F) = (l^,Curl$ w F). (105) 
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Here we use the notations 

l l = = e mIawM = eJawiA + e^dwlv] + e 3^M> ( 106 ) 

such that 

(i^, F) = r dw [x}F x + r dw [y}F y + r dw [z\F z . (107) 

The integral (11061) can be considered as a fractional line integral. 

In the right-hand side of fl!05j> . Ijy is a fractional surface integral over S = W such that 

I a s = I^ = e 1 I% r [y,z]+e a I%[z,x] + e 3 I% r [x,y]. (108) 

The fractional curl operation is 

CurlfyF = e l£lmn c D^[x m ]F n = ei ( c D^[y}F z - c D^[z\F y ) + 

+ e 2 ( c D^[z]F x - c D^[x]F z ) + e 3 ( c D% r [x]F v - C D^[y]F x ) . (109) 
For a = 1, equation (11091) gives the well-known expression 

Curl^F = curlF = eiEimn D Xm F n = e x (D y F z - D z F y ) + 

+ e 2 (D Z F X - D X F Z ) + e 3 (D x F y - D y F x ) . (110) 
The right-hand side of Eq. H105[) means 

(i^CWZ^F) = ( c D a w \y\F z - c D^[«]F y ) + 

+ I^[z,x] ( ^[z]^ - c ^[x]F z ) + I^[x,y] ( c Z^[x]F, - c D^[y]F x ) . (Ill) 
This integral can be considered as a fractional surface integral. 

8 Fractional Gauss's formula 

Let us give the basic theorem regarding the Gauss's formula in a fractional case. 
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Theorem (Fractional Gauss's Theorem for a Parallelepiped) 

Let F x (x,y,z), F y (x,y), F z (x,y,z) be continuously differentiable real-valued functions in a do- 
main that includes the parallelepiped 

W:={(x,y,z): a < x < b, c<y<d, g<z<h}. (112) 

If the boundary ofW be a closed surface dW, then 

(l% w ,F)=I%Div%F. (113) 

This equation can be called the fractional Gauss 's formula. 

Proof. 

For Cartesian coordinates, we have the vector field F = F x ei + F y e 2 + F z e 3 , and 

Iw = J w[ x > V' z l l dw = eitawlyi A + e 2 iawr[a:, A + ^awl x ^ v\- ( 114 ) 

Then 

(l^,F) =I% w \y,z]F 9 + I% w [x,z\F y + I% w [x,y}F z , (115) 

and 

I&Div^F = I^[x,y,z]( c D% w [x}F x + c D^ w [y]F y + c D a dw [z]F z ). (116) 
If W is the parallelepiped 

W:={a<x<b, c < y < d, g < z < h}, (117) 

then the integrals (11141) are 

Ifr[x, y, z] = X[x] J2[y] g I£[z], (118) 

and 

tt^cffl^], (119) 

tt^afflAl^ (120) 
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q w [z,y]= a I?[x] JM- (121) 
As a result, we can realize the following transformations 



I a aw , F j = I% w [y, «]F, + J^[*, x)F v + I% w [x, y]F z = 
My} 9 Ih[A {F x (b,y,z)-F x (a,y,z)} + a I?[x] g I%[z] {F y (x,d, z) - F y (x,c, z)} + 



+ a mJM [F z (x, y, g) - F z (x, y, h) 
= X[x] c F d [y] g I£[z} { c a D a x [x'}F x {x\ y, z) + c c ' D a y [y'}F y {x , y' , z) + ° D«[z']F z (x, y, z') 

= /^( c D^,f) =I«DivfrF. 
This ends the proof of the fractional Gauss's formula for parallelepiped region. 



Remark. To define the triple integral and the theorem for non-parallelepiped regions R, 
we consider the function f(x,y,z), that is defined in the parallelepiped region W such that 
R C W, such that 

\F{x,y,z), (x,y,z)eR; 
f(x,y,z) = \ (122) 

[o, (x,y,z)eW/R 

Then we have 

I%[x,y,z]F(x,y,z) = Ify[x, y, z]f(x, y, z). (123) 

As a result, we define a fractional triple integral over the non-parallelepiped region R, through 
the fractional triple integral over the parallelepiped region W. 

9 Fractional differential forms 

9.1 Brief description of different approaches 

A fractional generalization of differential has been presented by Ben Adda in [T6l [T7] . A 
fractional generalization of the differential forms has been suggested by Cottrill- Shepherd and 
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Naber in [22] (see also [23j El])- The application of fractional differential forms to dynamical 
systems are considered in [26], [27]. Fractional integral theorems are not considered. 

(1) In the papers [T6l [T7] , the fractional differential for analytic functions is defined by 



I 

T(l + a) 



d °f = W1 , ^ dxjN°f(x), (124) 



where N£ are Nishimoto fractional derivatives [8] (see also Section 22 of [2]), which is a gener- 
alization of the Cauchy's differentiation formula. 

(2) In the paper [22] (see also [251 I21H [21]). an exterior fractional differential is defined 
through the Riemann-Liouville derivatives by 

n 

d a = Y,(dx,r D a X] . (125) 

3=1 

In two dimensions n = 2, 

d a =(dx) a D« + (dy) a D«, 

suhc that 

ft = (i) "iR + ( *' , i^i' (126) 



where we use 

o 



1 = — J r. (128) 



(3) In the paper [27], an exterior fractional differential is defined through the fractional 
Caputo derivatives in the form 

n 

d a = J2(dx,r%D a Xj . (129) 

3=1 

For two dimensions (x,y), we have 



<r=(dx)"%D2 + (dy) a %D* 

such that 



x i- a 



d a x = (dxT— (130) 

1 (2 — a) 
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d a y = (dy) a -£ - (i3i; 

1 (2 — a) 



(compare with (I126p and (I127p ). Equation (I129p can be rewritten as 



d* = J2T(2-a)x';- 1 d a x^D a Xj . (132) 

3=1 



This relation is used in 1271 as a fractional exterior differential. 



9.2 Definition of a fractional exterior differential 

A definition of fractional differential forms must be correlated with a possible generalization of 
the fractional integration of differential forms. To derive fractional analogs of differential forms 
and its integrals, we consider a simplest case that is an exact 1-form on the interval L = [a, b}. 
It allows us to use the fractional Newton-Leibniz formula. 

In order to define an integration of fractional differential forms, we can use the fractional 
Riemann-Liouville integrals. Then a fractional exterior derivative must be defined through the 
Caputo fractional derivative. 

Equation (|67|) of FTFC means that 
r x dr' 

ifM^^' 1/1/(/) = /W " /(fl) ' (0 <«<!). (133) 

Using 

dx' = (dx') l ' a (dx') a , (0 < a < 1), 
equation fl 133[) can be presented in the form 

[ r(aKx-xr-° {( dx Ta Da AxV(x")) = fix) - /(a), (0 < a < 1). (134) 
The expression in the big brackets of fll34p can be considered as a fractional differential of the 
function f(x). As a result, we have 

I a L [x] a d a J(x) = fib) - /(a), (0 < a < 1), (135) 

where L = [a, b], and the fractional integration for differential forms is defined by the operator 

^ : =/idS^- (136) 
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The exact fractional differential 0-form is a fractional differential of the function 

a d%f{x):={dxY c a D a x [x']f{x'). (137) 

Equation (11351) can be considered as a fractional generalization of the integral for differential 
1-form. 

As a result, the fractional exterior derivative is defined as 

a d a x :={dx m r c a D a Xm [x' m }. (138) 
Then the fractional differential 1-form is 

u(a) = {dx m ) a F m {x). (139) 
The exterior derivative of this form gives 

a d%u{a) = (dx m ) a A (dx n ) a ^D: n [x']F m (x'). (140) 
To prove the proposition (11401) . we use the rule 

oo 

D a x {fg) = Y.^)^ D T s f)D s x g, 



s=0 



and the relation [5] 
for integer s, where 



D s [x](dx) a = (s>l), 
{k) r(l - a)F(k + 1) ' 



For example, we have 

oo 

d a i(dx m ) a F rn ] = J2(dxnY A ( a k ) (aD:; s [x' n }F m (x'))D s [x n ](dx m ) a 



s=0 



= {dx n ) a A (dx m ) a (JJ) D« n [x' n ]F m {x') = (D« n [x' n ]F m {x')) (dx n ) a A (dx m ) a . 
Using the equation (see Property 2.16 in [5]) 

a D«[x'](x' - af = f!f + 1) (x - a)"-*, (141) 
1 (p + 1 — a) 
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where n — 1 < a < n, and ft > n, and 

aD^fs'Kx' - a) k = (Jfe = 0, 1, 2, n - 1), 



(142) 



we obtain 



a d a x {x-a) a = (dx) a C a D«[x']x' 



(dx) a T(a + l). 



(143) 



Then 



and the fractional exterior derivative (11381) is presented 




(144) 




(145) 



The fractional differential 1-form (11391) can be written as 



u>(a) 



1 




(146) 



r(a + i) 



Remark. Using the suggested definition of fractional integrals and differential forms, it is 
possible to define a fractional integration of ra-form over the hypercube [0, l] n . Unfortunately, 
a generalization of this fractional integral, which uses the mapping of the region W C M n 
into [0, l] n , has a problem. For the integer case, we use the equation 



For the fractional case, the chain rule for differentiation (the fractional derivative of composite 
functions) is more complicated (see Section 2.7.3. [1]). As a result, a consistent definition of 
fractional integration of differential form for arbitrary manifolds is an open question. 

9.3 Differential vector operations through the differential forms 

To define a fractional divergence of the field F, we can consider the 2-form 



Dlfi^x)) = (Dlf)(Dlcf>). 



(147) 



uj 2 = F z dx A dy + F y dz A dx + F x dy A dz. 



(148) 
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Then the fractional exterior derivative of this form is 

duj 2 = (D X F X + D y F y + D z F z )dx A dy A dz = divF dx A dy A dz. (149) 
To define a fractional generalization of the curl operation for F, we can use the 1-form 

ui = F x dx + Fydy + F z dz. (150) 
Then the fractional exterior derivative of this 1-form is 

duj x = (D x F y - D y F x )dx A dy + (D y F z - D z F y )dy A dz + (D X F Z - D z F x )dx A dz. (151) 
To define the fractional gradient, we consider the 0-form 

Wo = f(x,y,z) (152) 
Then the fractional exterior derivative of / gives 

3 

duo = D x fdx + D y fdy + D z fdz = ^^(gradf)kdx k . (153) 

k=l 

It is not hard to obtain fractional generalizations of these definitions. 

10 Fractional nonlocal Maxwell's equations 
10.1 Local Maxwell's equations 

The behavior of electric fields (E, D), magnetic fields (B,H), charge density (p(t,r)), and 
current density (j(t,r)) is described by the Maxwell's equations 

divB(t,r)=p(t,r), (154) 

curlE(t,r) = -dtB(t,r), (155) 

divB(t,r) = 0, (156) 

curl H(t, r) = j(t, r) + d t T)(t, r). (157) 
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Here r = (x,y,z) is a point of the domain W. The densities p(t,r) and }(t,r) describe an 
external sources. We assume that the external sources of electromagnetic field are given. 
The relations between electric fields (E, D) for the medium can be realized by 

-+oo 



e(r, r')E(t,r')dr' 

-oo 



(158) 



where Eo is the permittivity of free space. Homogeneity in space gives e(r,r') = e(r — r'). 
Equation (j!58p means that the displacement D is a convolution of the electric field E at other 
space points. A local case corresponds to the Dirac delta-function permittivity e(r) = e5(r). 
Then Eq. flTgHJ) gives D(t,r) = e eE(t,r). 

Analogously, we have nonlocal equation for the magnetic fields (B,H). 

10.2 Caputo derivative in electrodynamics 

Let us demonstrate a possible way of appearance of the Caputo derivative in the classical 
electrodynamics. If we have 

/+oo 
e(x-x')E(t,x')dx', (159) 
-oo 



then 

-+oo 



/+oo r+oo 
[D l x e{x - a;')] E (^ x'W = ~ I \D l x ,e{x - x')]E(t, x')dx' . (160) 
-oo J — oo 

Using the integration by parts, we get 

/+oo 
e(x-x')Dl,E(t,x')dx'. (161) 
-oo 

Consider the kernel e(x — x') of integral (11611) in the interval (0,x) such that 

{e(x — x'), < x' < x; 
(162) 
0, x' > x, x' < 0, 

with the power-like function 

e(x-x')= 1 1 (0<a<l). (163) 

T(l — a) [x — x') a 

32 



Then Eq. ( 11611) gives the relation 

DlD(t, x) = oD%E(t, x), (0 < a < 1) (164) 
with the Caputo fractional derivatives q D%. 

10.3 Fractional nonlocal Maxwell's equations 

Fractional nonlocal differential Maxwell's equations have the form 

Div%E(t,r)=g lP (t,r), (165) 

Curl%E(t,r) = -d t B(t,r), (166) 
Div%B(t,r) = 0, (167) 
g 2 CurQ B(t, r) = j(t, r) + g^d t E(t, r), (168) 

where a s , (s = 1, 2, 3, 4), can be integer or fractional. 

Fractional integral Maxwell's equations, which use integrals of noninteger orders, have been 
suggested in [25] to describe fractional distributions of electric charges and currents. 

In the general form, the fractional integral Maxwell's equations can be presented in the form 



(l^,E(t,r)) = 9l Qp(t,r), (169) 

(l-,E(t,r))=-|(l-,B(t,r)), (170) 

(l^,B(t,r)) =0, (171) 

g 2 (l^ s ,B(t,r)) = (l? 4 ,j(t,r)) +^ 1 |(l^,E(t,r)). (172) 



These fractional differential and integral equations can be used to describes an electromag- 
netic field of media that demonstrate fractional nonlocal properties. The suggested equations 
can be considered as a special case of nonlocal electrodynamics (see j3TJ [32J, [331 Ell 135]). 

Fractional coordinate derivatives are connected with nonlocal properties of the media. For 
example, a power-law long-range interaction in the 3-dimensional lattice in the continuous limit 
can give a fractional equation [40J . 
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10.4 Fractional conservation law for electric charge 

Let us derive a conservation law equation for density of electric charge in the region W from 
the fractional nonlocal Maxwell's equations. 
The time derivative of (I169p is 

Div%a t -E(t,r)=g 1 d t p(t,r). (173) 

Substitution of (11681) into ( 1173[) gives 

9*Div% (g 2 Curl%B{t,r) - j(f,r)) = gi d t p(t,r). (174) 

If ai = Q!4, then 

£>ii$ CwZ$ B(t, r) = 0, (175) 

and we have the law 

gidtp(jk,r)+g 3 Div%j(t,r) = 0. (176) 

This fractional equation is a differential form of charge conservation law for fractional nonlocal 
electrodynamics. 

If ai = Q!4, we can define the fractional integral characteristics such as 

Qw(t) = gil%[x,y,z]p(t,x,y,z), (177) 

which can be called the total fractional nonlocal electric charge, and 

Jsw(t) = g 3 (j%y,i) =g 3 {l^[y,z\j x + I^[z,x]j y + I^[x,y\j z ) (178) 

is a fractional nonlocal current. Then the fractional nonlocal conservation law is 

j t Qw(t) + J dw (t) = 0. (179) 

This integral equation describes the conservation of the electric charge in the nonlocal electro- 
dynamics for the case a% — 014. 
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10.5 Fractional waves 

Let us derive wave equations for electric and magnetic fields in a region W from the fractional 
nonlocal Maxwell's equations with j = and p = 0. 
The time derivative of Eq. (11661) is 

d 2 B = -CurQd t E (180) 

Substitution of fTT72D and j = into (1T80]) gives 

«9 2 B = -g 2 g 3 Curl% CurQB{t,r). (181) 

Using flHDl and (11671) for ct 2 = a 3 = a^, we get 

«9 2 B = ^ 3 ( C D^) 2 B. (182) 

As a result, we obtain 

<9 2 B-w 2 ( c D^) 2 B = 0, (183) 

where v 2 = g 2 g3- This is the fractional wave equation for the magnetic field B. Analogously, 
Eqs. (I166P and (I172p give the fractional wave equation for electric field 

d t 2 E--u 2 ( c D^) 2 E = 0. (184) 

The solution B(t,r) of equation H183[) is a linear combination of the solutions B + (t,r) and 
B_(t,r) of the equations 

d t B + (t,r)-v c '-D^B + (t,r) = 0, (185) 

<9 t B_(t,r) + v c T>^B_{t,r) = 0. (186) 

As a result, we get the fractional extension of D'Alembert expression that is considered in [36J. 
For the boundary conditions 

lim B(t,r) = 0, B(t,0) = G(t), (187) 

the general solution of equations fl 185[) and (11861) is given [5] by 

i r+ca 

B m± (t,r) = — / duE ail [Tivux^\ G m {w)e**, (188) 
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where G m {uS) = !F[G m (t)], and E a ^[z\ is the biparametric Mittag-LefHer function [5]. Here 
B± m (t,r), and G m (t) are components of B±(i, r) and G(t). 

For one-dimensional case, B x (x,y, z,t) = u(x,t), B y = B z = 0, and we can consider the 
fractional partial differential equation 

D 2 u(x, t) - v 2 D 2 x a u{x, t) = 0, x G E, x > 0, v > 0, (189) 

with the conditions 

D k x u(0,t) = f k (t), (190) 

where k = for < a < 1/2, and A; = 1 for 1/2 < a < 1. If < 2a < 2 and u > 0, the system 
of equations (11891) . (11901) is solvable (Theorem 6.3. of [5]), and the solution u(x, t) is given by 

u(x,t) = J2 G 2 k a (y,t)f k (y)dy, (n-Ka<n), (191) 

fc=0 

where 

G 2 k a (x, t) = ^vx k ~ a (p(-a, k + l-a, v\t\x~ a ). (192) 

Here 0(— a, k + 1 — a,t> |t|x _a ) is the Wright function [5]. 

Note that the solutions of equations as fl!85[) and (11861) are based primary on the use 
of Laplace transforms for equations with the Caputo q D™ derivatives. This leaves certain 
problems [5] with the fractional derivatives ^D™ for a G R. 

11 Conclusion 

Let us note some possible extensions of the fractional vector calculus. 

(1) It is very important to prove the suggested fractional integral theorems for a general 
form of domains and boundaries. 

(2) It is interesting to generalize the formulations of fractional integral theorems for a > 1. 

(3) A proof of fractional theorems for differential forms can be interesting to formulate a 
fractional generalization of differential geometry. 
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In the fundamental theorem of fractional calculus (FTFC) we use the Riemann-Liouville 
integration and the Caputo differentiation. The main property is that the Caputo fractional 
derivative provides us an operation inverse to the Riemann-Liouville fractional integration from 
the left. Note that a fractional generalization of the differential vector operations and the 
integral theorems for the fractional integro-differentiation of Riesz, Grunvald-Letnikov, Weyl, 
Nishimoto is an open problem. 

There are the following possible applications of the fractional variational calculus (FVC). 

(a) A fractional nonlocal electrodynamics that is characterized by the power law non-locality 
can be formulated by using the FVC. 

(b) Nonlocal properties in classical dynamics can be described by the FVC and by possible 
fractional generalizations of symplectic geometry and Poisson algebra. In general, fractional 
differential forms and fractional integral theorems for these forms can be used to describe 
classical dynamics. 

(c) A possible dynamics of fractional gradient and Hamiltonian dynamical systems can be 
described by the FVC. 

(d) The continuum mechanics of fluids and solids with nonlocal properties (with a nonlocal 
interaction of medium particles) can be described by the FVC. 

The fractional derivatives in equations can be connected with a long-range power-law inter- 
action of the systems [37J [38], HQ] . The nonlocal properties of electrodynamics can be considered 
[39J as a result of dipole-dipole interactions with a fractional power-law screening that is con- 
nected with the integro-differentiation of non-integer order. For noninteger derivatives with 
respect to coordinates, we have the power-like tails as the important property of the solutions 
of the fractional equations. 
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